Transport coefficients of soft sphere fluids at high densities 
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Molecular dynamics computer simulation has been used to compute the self-diffusion coefficient, 
and shear viscosity of soft-sphere ffuids, in which the particles interact through the soft-sphere or 
inverse power pair potential. The calculations have been made along the melting line in a wide range 
of pressures and temperatures. The validity of scaling relations for thermodynamic parameters and 
kinetic coefficients was checked. It was shown that the Stokes-Einstein relationship is obeyed if the 
Barker diameter is used as a characteristic length scale. It was also shown that the viscosity is 
non-monotonic along the isochores as predicted by Ya. Rosenfeld. It was shown that the viscosity is 
strongly growing along the melting line, however, this increase does not stimulate the glass transition 
because the relaxation time is decreasing. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 
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I. INTRODUCTION 

Knowledge of transport coefficients is extremely im- 
portant for characterizing the properties of liquid. They 
allow one to characterize how fast the dynamics of the 
liquid is. Because of this, transport coefficients are inten- 
sively studied experimentally, theoretically and by sim- 
ulations. The most important transport property for an 
experimental study is viscosity. Viscosity can easily be 
measured by a number of techniques. Nevertheless many 
questions on the behavior of viscosity as a function of 
pressure and temperature still remain. Some of them are 
considered in the present article. 

It is well known that the viscosity of most liquids grows 
exponentially with increasing pressure at constant tem- 
perature and demonstrates critical behavior in the vicin- 
ity of the glass transition line. In the case of hard or 
soft spheres the glass transition at high pressures cor- 
responds to the jamming of the system. However, such 
kind of behavior is typical for supercooled liquids while 
the behavior of viscosity close to the melting line is not 
clear. 

As was discussed in [1], there are two empirical ap- 
proaches to the description of the viscosity behavior of 
simple liquids under simultaneous changes of tempera- 
ture and pressure. The first approach goes back to Bridg- 
man's works ^ and suggests that the viscosity of the 
liquid is nearly constant along the isochores. The other 
and currently more widespread approach was formulated 
by Fourier [3|. It is based on the presumption that 
the viscosity of the melts is invariable along the melting 
curve. The experimental data analysis shows that both 
approaches are incorrect; the constant viscosity lines for 
rare gas liquids, as well as for liquid metals, have a slope 
that is intermediate between the melting line and the iso- 
chore slopes jl (figures). As a consequence, the viscosity 
of simple melts along the melting curve increases and this 
increase can be extrapolated to megabar pressures [1]. 
However, this extrapolation is in fact difficult because of 



a rather small viscosity variation in the covered range [HI . 
For example, the viscosity of the Fe melt grows along the 
melting curve by several times with the pressure increase 
to 10 — 15 GPa [1, 0]- However, the obtained data do 
not permit drawing any conclusions with regard to the 
character of this growth; various analytic models equally 
closely descriptive of the studied baric dependencies give 
different extrapolation results. In consequence, different 
equations employed in the extrapolation of the Fe melt 
viscosity to pressures of 1.4 — 3.1 Mbar, corresponding to 
the conditions existing in the Earth's outer core, provide 
the data that vary over 10 — 20 orders of magnitude Q 
(figure 6). The same is valid for other simple melts like 
liquid Ar. 

Thus, it is presently impossible to draw inferences 
about how high the viscosity can grow along the melt- 
ing curve for the given melts with the pressure increase 
to about 100—1000 GPa. If the viscosity growth is minor, 
liquid metals and rare gas liquids will remain very poor 
glass formers at megabar pressures. On the other hand, 
the viscosity growth along the melting curve amounts 
to several orders of magnitude, so the question arises 
whether the liquid metals and rare gas liquids become 
viscous glass-forming systems under the increasing pres- 
sure. In this case, the glass transition of the given melts 
in the megabar range should be similar to the jamming 
of soft spheres. 

The main purpose of this paper is the computer simu- 
lation study of the behavior of the transport coefficients 
of a soft sphere system in a wide range of densities and 
temperatures. Although the interatomic potential in liq- 
uid metals can not be described by a pair potential, the 
potential in rare gas liquids can be approximated by the 
Lennard-Jones function which can be substituted in the 
high density and temperature limit by a soft sphere po- 
tential. In this respect it is interesting to carry out a sys- 
tematic study of the transport coefficients of soft sphere 
system as a generic model for real liquids at high pres- 
sure. 
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Note that most of previous studies of transport coeffi- 
cients of liquids were concerned to the low density - low 
temperature region while high density - high tempera- 
ture properties were beyond the consideration. Only few 
works on the high density transport properties are avail- 
able (see, for example, [7|) In this sense the present work 
is directed to fill this gap. 



This approach was further developed in the work [12 



where the scaling relations were derived under the as- 
sumption of scaling invariancy of the Hamilton equations 
for homogeneous potentials. 

Remember that kinetic energy is a homogeneous func- 
tion of the power 2: 



II. SOFT SPHERES: THEORETICAL APPROACH 

Interestingly, one can make some important predic- 
tions about thermodynamic properties of soft spheres on 
the basis of the form of the potential only. To do this 
one can use the Klein theorem [9, 10]. 
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Consider the dimensionless density and reduced 
coordinates s = r{^y/^. The partition function of soft 
spheres in these variables has the following form [11|]: 



(6) 



If scaling = ar^, pi = hpi 



and 
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is applied then the new phase trajectory of 



the system has the same geometric shape [l^. Then the 
following scaling laws are valid (we skip the Boltzmann 
constant in the formulas below for the sake of brevity) : 



V = a-'Va 



(7) 
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where A = — -, — is de Broile wave length and 

13 = l/{kBT). The Klein theorem states that if t/ is a 
homogeneous function of degree n, i.e. U (Ar i , Arjv) = 
A"[/(ri, rjv) then for N sufficiently large [lO| 



ea:p[-/3-[/(ri,...,rjv)]rfri...drjv = G(/3-p-"/3)^ 

(2) 

where G is some function of (/3p "Z^). 

Using the Klein theorem one can find that the equation 
of state has the form: 
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Using the Green - Kubo relations for the diffusion 
D = /q°° < V(t)V(0) > dt and viscosity 77 = lo" < 

a'^'^ {t)a^y (0) >, where a^^ is a stress tensor component, 
one can derive the scaling relations for diffusion and vis- 
cosity: 



P2n + e 



(10) 



or 



PV 
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= l + 0(p(^f/") 



(3) 



If the system with such equation of state has a first 
order phase transition (let us call it melting), then the 
transition line is characterized by the universal parameter 
p(jpy)^/" 11 1 . From this it follows that the densities of 



a liquid and crystal along the melting line are: 



and 



or 



77 



^ p 2(^ + 3) 
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where ci and Cc are constants. One can see that the 
densities along the melting line follow some scaling law. 



Below we compute different quantities in simulation 
and fit the results to the corresponding scaling law using 
the least square procedure. 



One can see that, following [1^, the scaling relations 
can be generalized to be applicable to the transport co- 
efficients, however this approach involves some assump- 
tions which are not fully proved. This is why it is particu- 
larly interesting to check whether the scaling predictions 
for the transport coefficients are correct. 
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FIG. 1: Pressure as a function of temperature along the melt- 
ing line in usual (a) and scaled (b) forms. The symbols are 
MD data and the line is the scaling law fitting. 



III. SIMULATION 

In the present work we study the viscosity, rjs , and dif- 
fusivity, D, of soft spheres. The choice of the system is 
determined by the following reasons. In order to investi- 
gate the behavior of the system along the melting line in 
a wide range of pressures we need to know the phase di- 
agram of the system including a high pressure region. In 
most cases the melting lines were studied in the vicinity 
of the triple points. As it was shown that the soft sphere 
system obeys the scaling law. The parameters of the scal- 
ing of pressure and density with temperature along the 
melting line were determined both from simulation 13 1 
and from density functional theory [14{ . 

The soft sphere potential has the form: 



$(r) 



(14) 



In our simulations we use n = 12. The potential is 
cut at the distance Tc = 2.5 as in standard simulations of 
Lennard-Jonse systems. Reduced units (e = 1 and u = 1) 
are used in the paper. As it was shown (see, for example, 
previous section and the behavior of soft spheres 

system depends on the parameter 7„ = p(T'^(^^)^^/". 
In 13| it was shown that the line of fluid to solid tran- 
sition is 7„ = 1.15. From this formula one can compute 



the transition density to be used in simulations. 

The system of 1000 particles was simulated in the mi- 
crocanonical (constant N, V and E) ensemble. Equations 
of motion were integrated by the velocity Verlet algo- 
rithm. The equilibration time was set to 5 • 10^ time 
steps and the production time - 3.5 • 10^ steps. The time 
step we use is dt — 0.0005. During equilibration the 
velocity rescaling was applied to keep the temperature 
constant. During the production cycle, the NVE en- 
semble was used. The diffusion coefficient was calculated 
from the mean square displacement via Einstein relation, 
while the viscosity was determined by the integration of 
the stress correlation function [§]. 

The computation of the phase diagram of the system 
was also done. The reason for this computation is as fol- 
lows. To our knowledge, all previous works studying soft 
sphere phase diagram considered the system at temper- 
atures of the order of unity. Although it allows finding 
the scaling coefficient 7„, the value may contain some 
errors. Since in the present work we study the temper- 
atures up to T = 100, these small errors can give quite 
large uncertainty in the location of the melting line. 

Because of this we computed the melting densities at 
T = 20 and T = 50. The same system of 1000 parti- 
cles was used for this calculation in the case of the liquid 
and of 1372 particles for the FCC solid. The transition 
points were determined via the double tangent construc- 
tion to the free energy lines. The crystal free energy was 
computed by coupling to the Einstein crystal Isl, while 
the liquid free energy was determined by integrating the 
equation of state from the dilute gas limit [8|. Table 1 
presents a comparison of the melting points obtained by 
our simulations to the predictions based on the scaling 
factor from Ref. [ic 
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pi (this work) 


Pi (ref. [13]) 


relative discrepancy 


20 


2.45 


2.43 


0.8% 


50 


3.12 


3.06 


1.6% 



We can see that the agreement is rather good. It con- 
firms that one can use the scaling parameters from the 
reference [131. 



IV. RESULTS AND DISCUSSION 
Transport Coefficients Beliavior 

The main purpose of this paper is to study the be- 
havior of transport coefficients - shear viscosity and 
diffusion - in a wide range of densities and tempera- 
tures. To do this, we measure transport coefficients 
along the melting line and along a set of isochores 
for temperatures ranging from T = 0.1 up to T = 
100.0. The following isochores are considered: p = 
0.65; 0.95; 1.51; 2.43; 2.69; 2.89; 3.06; 3.2; 3.33; 3.44. Note 
that the density and temperature dependence of pressure 
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FIG. 2: Viscosity along the melting line as a function of pres- 
sure and temperature. The symbols are the MD data and the 
line is the scaling law fitting. 



is rather strong; thus, the investigated density and tem- 
perature range corresponds to the pressure range varying 
from about unity to about 7000. Changing Lennard- 
Jones reduced units to real ones, we have the pressure 
varying in the range up to SOOGPa. It covers almost 
the entire range of experimentally attainable static high 
pressures. 

Figs, la and lb present the pressure dependence along 
the melting line in usual and scaled forms. One can see 
that the scaling law works well. 

The next figures (Figs. 2a and 2b) show the behavior 
of viscosity along the melting line as functions of both 
pressure and temperature and the fitting to the corre- 
sponding scaling formulas. As one can see, the scaling 
rules for viscosity are well fulfilled too. Unlike some pre- 
dictions, the viscosity is not constant along the melting 
line. The mechanism of the viscosity growth will be dis- 
cussed below. 

Finally, Figs. 3a and 3b present the diffusion coefficient 
along the melting line. 

Figs. 2 and 3 show the diffusion and viscosity growth 
along the melting line. Although it may seem counter 
intuitive, they satisfy to the scaling laws. We will discuss 
this problem below. 
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FIG. 3: Diffusion coefficient along the melting line as a func- 
tion of pressure and temperature . The symbols are the MD 
data and the line is the scaling law fitting. 



Figs. 4a and 4b show the diffusion coefficient along dif- 
ferent isochores and along the melting line as functions 
of temperature and pressure. As one can expect the dif- 
fusion increases upon temperature growth while at fixed 
temperature it decreases with increasing density. Dif- 
fusion also increases with increasing pressure along the 
isochores because in this situation we have simultaneous 
increase of temperature. 

The behavior of viscosity along the isochores is more 
complex (Figs. 5a and 5b). One can see from the in- 
sets in the figures that the viscosity is not a monotonic 
function of pressure and temperature. The viscosity de- 
creases with pressure close to the melting line, and in- 
creases away from it. Another conclusion from these fig- 
ures is that the viscosity behaves qualitatively similar 
along different isochores - the shapes of the curves are 
almost the same for different densities. 

A possible explanation of the nonmonotonic behavior 
of the viscosity was suggested by Rosenfeld In his 
works ^Td, li'] Rosenfeld showed the connection between 
transport coefficients of a fluid and excess entropy. Con- 
sidering the behavior of the entropy he showed that the 
viscosity can have a minimum. He also gave a qualita- 
tive explanation of this fact. If we consider viscosity as 
an integral over stress correlation function, we can split 
it to three contributions: kinetic - kinetic (fcfc), kinetic - 
potential (kp) and potential - potential fpp)|8l. Il7|: 
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I Ckkit)dt+ / Ckp{t)dt+ / Cpp{t)dt, (15) 
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where 
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where Vx,i{t) denotes x component of the velocity of the 
particle i at time t. 
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where Xjk{t) and yjk(t) are Cartesian components of 
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The role of different contributions to the viscosity for 
a Lennard - Jones fluid was studied by several authors 
(see, for example, (isl. K^). 

One can expect that close to the melting line the main 
contribution to the viscosity comes from the potential 
- potential correlation, while at high temperatures the 
kinetic - kinetic part should be dominant. Having the 
pp correlation decreasing and the kk one increasing, we 
can expect that a minimum of the total viscosity can be 
observed. 

In order to verify this suggestion we carried out the 
calculations of the different contributions to the viscos- 
ity for the density p = 1.51 and a set of temperatures 
(pressures) (Figs. 6a and 6b). As one can see from these 
figures the suggestion is correct: while the pp contribu- 
tion decreases with increasing temperature both the kk 
and kp ones increase. Because of this one can expect the 
further increase of the viscosity with increasing temper- 
ature. 

Given this increase of the viscosity, one can pose a 
question about a glass transition in the soft sphere sys- 
tem. Two definitions of the glass transition point are 
common: the one which is connected to the relaxation 
time of the system, and another which defines the glas- 
sification through the viscosity increase. The Maxwell 



relaxation time and shear viscosity are related via the 
equation [17[: 



where 



Ginf = pksT - 



15 



^^9ir)^^r -) 



(19) 



(20) 



is the infinite-frequency shear modulus lO, 21 1 . 

The behavior of the elastic modulus Gmf along the 
melting line is shown in Fig. 7. One can see that Ginf 
grows dramatically. Fig. 8 presents the relaxation time 
along the melting line. As it is clear from this plot, 
although the viscosity of the system increases with in- 
creasing temperature along the melting line the relax- 
ation time becomes very small. Having small relaxation 
time, the system can definitely not become glasslike. 

One can conclude that the main glassification criterion 
is the one through the relaxation time. As a result, even 
at certain conditions one has the viscosity increase, at the 
same time the high viscosity does not mean glassification 
because of decrease of the relaxation time. It seems that 
it is a high pp viscosity that brings the system to glass 
transition while dominating kk viscosity does not. 
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FIG. 4: Diffusion coefficient along several isochores and the 
melting line as functions of temperature and pressure. (Color 
online) 



Stokes - Einstein Relation 



FIG. 5: Viscosity along the isochores (squares) and along the 
melting line (stars)as a function of pressure and temperature. 
The symbols are MD data and the line is the scaling law fit- 
ting. The insets represent the low pressure or low temperature 
behavior of the plots. (Color online) 



Given the behavior of the viscosity and diffusion coef- 
ficient, it is interesting to consider the Stokes - Einstein 
relation. It can be written as 



TTaDrj ' 



(21) 



where a is the diameter of a particle and c is the constant 
which equal 2 for "stick" boundary conditions and 3 for 
the "slip" ones [l7|- Several authors studied the Stokes 
- Einstein coefficient c in different systems (22| - [26t . The 
deviation from the Stokes - Einstein law was found. In 
Reference ^3^, a modified relation was suggested. To 
take into account the effect of a change of the diameter 
a, it was suggested to use the relation: 



knT 



TraDrjp 



1/3- 



(22) 



Although this normalization reduced the deviation, it 
did not make the coefhcient c constant. The study of 
the Stokes - Einstein relation in the soft sphere system 



24 - 26i also confirmed the deviation from both above for- 



mulas. 

Here we apply the scaling laws to understand the 
Stokes - Einstein coefficient behavior. Recall that along 
the melting line the diffusion coefficient changes as Z? ~ 
T^/i^ while ry ~ T'^/^. Inserting it in Eq. ^ we ob- 
tain c ^ 2^-1/12^ Figs. 9a and 9b represent the Stokes - 
Einstein coefficient c along the melting line in usual and 
scaled forms. 

One can see that scaling law works for the Stokes - 
Einstein coefficient too [27| . The apparent choice for the 
particle size a is the Barker diameter. For soft sphere 
system the Barker diameter is defined as 
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FIG. 6: Different contributions to viscosity as a function 
of pressure and temperature along the isochore p = 1.51. 
Squares - full viscosity, circles - pp contribution, up triangles 
- kk contribution and down triangles - kp contribution. (Color 
online) 
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FIG. 7: Infinite- frequency shear modulus dnf along the melt- 
ing line. 
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(23) 



i.e. for our case c ~ y-1/12 Yvoto. this one can conclude 
that the Barker diameter should be used as a character- 
istic particle size. Fig. 10 shows c/T^/^^. One can see 
that the Stokes - Einstein coefficient becomes constant 
and equals 2 with this normalization. 



Isoviscosity Lines 

As it was mentioned in the introduction, the shape 
of the isoviscosity lines is still unclear. There are two 



FIG. 8: Relaxation time along the melting line. 



suggestions claiming that the viscosity is constant along 
the isochore or along the melting curve, but both of them 
contradict to the experimental data [J]. In this respect 
it is interesting to see the behavior of isoviscosity lines 
in the soft sphere system as a generic model for more 
complex liquids. 

Figs. 11 and 12 present the isoviscosity lines in com- 
parison with the viscosity along the melting line and 
along the isochores. One can see that the slope of the 
isoviscosity lines is higher then the viscosity along the 
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FIG. 9: Stokes - Einstein coefficient c in usual (a) and scaled 
(b) forms. Squares - MD results, line (b) is linear fit. 
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FIG. 11: Isoviscosity lines compared with the viscosity along 
the melting line and isochores. Squares - melting line, opened 
triangles - isochores (from left to right: p — 2.43, 3.06 and 
3.33), filled circles - isoviscosity lines (from left to right: i] = 
41.61, 55.32, 77.94 and 88.86). (Color online) 
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FIG. 10: Stokes - Einstein coefficient cjT^I^^ . 



melting line, but lower then the one along the isochores. 
One of us discussed the behavior of the isoviscosity lines 
in the previous publication [5l • Based on the experimen- 
tal data for mercury and argon, it was suggested that the 
isoviscosity lines are located between the melting line and 
isochore. From this one can see that the simulation re- 
sults qualitatively correspond to the experimental data. 



Excess Entropy Scaling 



FIG. 12: Isoviscosity lines compared with the viscosity along 
the melting line and isochores. Squares - melting line, 
filled circles - isoviscosity lines (from left to right: r\ = 
41.61,55.32,77.94 and 88.86). Note that the isochores are 
vertical lines in these coordinates (dashed lines here). (Color 
online) 



namic and thermodynamic properties of liquids. One of 
the most common relations was developed by Rosenfeld 
lil 16| . According to Rosenfeld, one can define some re- 
duced transport coefficients which exponentially depend 
on the excess entropy of the liquid. The reduced diffusion 
D* and viscosity r\* have the following forms: 



It is important from both theoretical and practical 
points of view to establish some relationships between dy- 



D* =D 



,1/3 



(fcBr/m)i/2 



(24) 
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FIG. 13: Full (squares) and pair (circles) excess entropy of the 
soft spheres along the isochore p = 1.51. The inset represents 



the relative error —. 



in percent. (Color online) 



FIG. 14: Logarithm of reduced diffusion D* as a function of 
full (squares) and pair (circles) excess entropy. (Color online) 
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and their dependence on the excess entropy Sex = {S 
Sid)/{NkB) looks like 



X — ax ■ e 



bxS, 



(26) 



where X is the transport coefficient (D or 77 in the present 
case), and ax and hx are the constants which depend on 
the studied property |15l. Il6|. 

Very often one can replace the excess entropy by a 
pair contribution to it, which can be computed via the 
relation 



S2 



^-\p I dr[g(r)ln(g(r))-(g(r)-l)], 



(27) 



where g{r) is radial distribution function of the liquid. 

In the present work we verify the scaling relations for 
the diffusion coefficient and viscosity. To do this, we 
calculate the pair and full excess entropies along the 
isochore p = 1.51. The pair excess entropy is calcu- 
lated from the radial distribution function. For com- 
puting the full excess entropy, the following procedure 
is applied. As we mentioned above, the free energies of 
the system at T = 20.0 were calculated for determin- 
ing the melting point. After that we computed the en- 
ergy along the isochore p = 1.51. Having free energy at 
Tq = 20.0 and the temperature dependence of internal 
energy, one can find the free energy along this isochore: 
1^ ~ fclk = -^To d{ju)U{T'). The excess entropy along 
the isochore is obtained by differentiating the free energy 
with respect to temperature. 

Fig. 13 presents the excess entropy and pair contribu- 
tion to it along the isochore p = 1.51. One can see that 
at some temperature the relative difference of the en- 
tropies becomes approximately constant. In the case of 
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FIG. 15: Logarithm of reduced viscosity 77* as a function of 
full (squares) and pair (circles) excess entropy. (Color online) 



the exponential dependence, it gives a constant shift of 
the curves. 

Figs. 14 and 15 show the logarithm of the Rosenfeld 
scaled diffusion and viscosity along this isochore. As can 
be seen from these plots, the scaling relation is fulfilled 
with good accuracy in a wide range of temperatures. The 
curves corresponding to the full excess entropy and to the 
pair contribution to it have the same shape but shifted 
with respect to each other, which corresponds to the in- 
tuitive conclusion described above. 

In conclusion to this subsection, we would like to make 
some remarks. First of all, recently another mechanism 
for the orig in of the excess entropy scaling has been sug- 
gested [30|. According to this work, the scaling law 
is valid for the liquids with strong energy-virial corre- 
lations (the authors call them "strongly correlated liq- 
uids"). The class of strongly correlated liquids includes 
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many different systems, like, for example, a Lennard- 
Jones liquid, Kob- Andersen mixture, square-well liquid 
and many others. The soft spheres considered in this 
article are exactly the strong correlated liquids [31 1. 

On the other hand, the applicability of the entropy 
scaling to not-strongly-correlated liquids is questionable. 
It seems that in the case of the existence of thermody- 
namic anomalies in the system 3^, [s^ , Rosenfeld scal- 
ing does not work. In Ref. |3J it was shown that for 
Herzian spheres. Gauss Core Model, and soft repulsive 
shoulder potential, the Rosenfeld entropy scaling breaks 
down. These systems demonstrate the diffusion anoma- 
lies at low temperatures: the diffusion increases with in- 
creasing density or pressure. It is shown that for the 
first two systems, which belong to the class of bounded 
potentials, the Rosenfeld scaling formula is valid only in 
the infinite temperature limit where there are no anoma- 
lies. For the soft repulsive shoulder, the scaling formula 
is valid already at sufficiently low temperatures, however, 
out of the anomaly range. 



Influence of Attraction 
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Finally we would like to study the influence of attrac- 
tive forces on the behavior of the viscosity. To do this, 
we simulate the Lennard- Jones like system. The system 
has the potential function: 

<^>,j = s-[{^r-{^r] (28) 
r r 

with the parameters s — I and (7 = 1. 

We measure the viscosity of both LJ systems at the 
isochore p = 1.51. The resulting curves are shown in Figs. 
16a and 16b. One can see from these figures that the 
difference takes place only in the low temperature region 
and should be concerned to the pp contribution while at 
high temperatures where kk contribution is dominant the 
viscosity of soft spheres and LJ system have very close 
magnitudes. 

Based on these results, one can expect that in the case 
of real liquids the viscosity can demonstrate complex be- 
havior on compression, since they are characterized by 
complex non pair potentials which can even be state de- 
pendent. The results for the soft sphere system can di- 
rectly be applied only for the rare gas liquids, although 
even in this case the attraction is important, as is seen 
from Figs. 16a and b. 

As one can see from the figures the difference between 
the soft spheres and Lennard- Jones system takes place 
only in the low temperature regime. One can attribute 
it to the difference in potential - potential contribution 
to the viscosity while at high temperature the influence 
of the kinetic - kinetic term is not so strongly affected by 
the interatomic potential. 



FIG. 16: Viscosity of the soft sphere and Lenard- Jones sys- 
tems along the isochore p = 1.51. Squares - soft spheres, and 
circles - LJ. (Color online) 

IV. CONCLUSIONS 

In conclusion we have carried out a comprehensive 
study of the viscosity and diffusivity of the soft sphere 
system. We showed the validity of the scaling laws along 
the melting line. We proposed the scaling law for the 
Stokes - Einstein coefficient and verified it. The influence 
of different contributions to the viscosity was discussed. 
It was shown that the viscosity is strongly growing along 
the melting line, however, this growth does not stimu- 
late the glass transition because the relaxation time is 
decreasing. We calculated the shape of the isoviscosity 
lines and compared them with the viscosity along the 
melting line and isochores. The results we obtained are 
qualitatively similar to the experimental findings. The 
validity of the Rosenfeld excess entropy scaling was also 
verified. Interestingly, although the qualitative behavior 
of viscosity and diffusivity is rather different, the excess 
entropy scaling describes both coefficients with high ac- 
curacy. Finally, the influence of the attractive forces was 
verified via the comparison to the Lenard- Jones system. 
It was shown that the attraction is important only at 
sufficiently low temperatures. 
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